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Abstract 

In this paper, we study the Cauchy problem of a time-dependent drift-diffusion- 
Poisson system for semiconductors. Existence and uniqueness of global weak solutions 
are proven for the system with a higher-order nonlinear recombination-generation rate 
R. We also show that the global weak solution will converge to a unique equilibrium as 
time tends to infinity. 
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We consider the following drift-diffusion-Poisson model for semiconductors that is a cou- 
pled system of parabolic-elliptic equations: 
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n t = div(Vn + nV(ip + V^)) — R(n,p, x), 

P t = dw(Vp + pV(-ip + V p )) -R(n,p,x), (1.1) 
—e 2 Aip = n — p — D(x). 

System (jl.ip models the transport of the electrons and holes in semiconductor and plasma 
devices (cf. [I3[T8]). n = n(x,t) is the spatial distribution of electrons (negatively charged) 
and p = p(x,t) is the spatial distribution of holes (positively charged), ip = ifj(x,t) is 
the self-consistent electrostatic potential created by the two charge carrier species (electrons 
and holes) and by the doping profile D = D(x) of the semiconductor device. The charge 
carriers are assumed to be confined by the external potentials V n and V p . This replaces 
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the usual assumption of a bounded domain (cf. [HIIIHIIIB] and the references therein). The 
function R = R(n,p,x) represents the so-called recombination-generation rate for electron 
and hole. The parameter e appearing in the Poisson equation is the scaled Debye length of 
the semiconductor device that stands for the screening of the hole and electron particles. In 
this paper, we are interested in the Cauchy problem to system (II. ip and assume that (II. ip is 
subject to the following initial data 

n(x,t)\ t=0 = n I (x) >0, p(x,t)\ t=0 = pi(x) > 0. (1.2) 

The generation and recombination of electrons and holes in a semiconductor play an 
important role in their electrical and optical behavior [T2]. Recombination is a process by 
which both carriers annihilate each other: the electrons fall in one or multiple steps into 
the empty state that is associated with the hole. Generation can be viewed as its inverse 
process whereby electrons and holes are created. There are several typical recombination 
mechanisms that the energy of carriers will be dissipated during these processes by different 
ways (cf. e.g., [T2" | [TT 1 [T8] ) . For instance, 

1. Band-to-band recombination (also referred to as direct thermal recombination). The 
energy is emitted in the form of a photon. The recombination rate depends on the 
density of available electrons and holes and it can be expressed as 

R{n,p) = C(np-nf), (1.3) 

where rtj denotes the intrinsic carrier density of the semiconductor. 

2. Shockley-Read-Hall (SRH) recombination (also called the trap-assisted recombination). 
A two-step transition of an electron from the conduction band to the valence band 
occurs and R is in the form 

R(n,p)= . (1.4) 

r\n + r 2 p + r 3 

where ri,V2,V3 are proper positive functions. 

3. Auger recombination. An electron and a hole recombine in a band-to-band transition, 
but the resulting energy is given off to another electron or hole in the form of kinetic 
energy. The corresponding recombination rate is similar to that of band-to-band re- 
combination, but involves a third particle: 

R(n,p) = (C n n + C p p)(np - nf). (1.5) 

Extensive mathematical study of the drift-diffusion-Poisson system has been developed 
in the literature. For the initial boundary value problem of (jl.ip in a bounded domain 

C R N with various boundary conditions (e.g., the Neumann type or no-flux boundary 
conditions), existence and uniqueness as well as long-time behavior have been investigated 
by many authors, see for instance, [3Tf5tl7ti9l [T3l - [T5| [26] and reference therein. 
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For the sake of modeling simplicity and for the particularly interesting mathematical fea- 
tures, it is also interesting to consider the Cauchy problem of (11. ip . Existence and uniqueness 
results and stability of strong solutions in LP(R. N ) spaces (N > 2) were proven in [20] for a 
system analogous to (jl.ip . However, in their system there were no external potentials and 
the recombination-generation rate R was replaced by a given function / = f(t,x) that ex- 
pressed the variation of the charge by the external current. As far as the long-time behavior 
of global solutions to the Cauchy problem is concerned, when the recombination-generation 
term R is absent, exponential convergence to equilibrium with a confining potential and an 
algebraic rate towards a self-similar state without confinement have been obtained in pQ (cf. 
also [1]). The analysis therein is based on the well-known entropy approach for diffusion 
and diffusion-convection equations that has been extensively studied in recent years (cf. [21[6] 
and the references therein). We also refer to [19] in which an optimal LP decay estimate 
of solutions was obtained via a time weighted LP energy method (without confinement and 
recombination-generation rate R). When the recombination-generation process is taken into 
account, the situation is more complicated. In [23], the authors proved the global existence 
and uniqueness of weak solutions of problem (jl.ip in IR 3 with an (unbounded) external po- 
tential V n = V p = V and under the restrictive assumption that R has a linear growth (which, 
however, recovers the Shockley-Read-Hall recombination, cf. (jl.4p ). Besides, existence and 
uniqueness of the steady state and partial result on the convergence to equilibrium were ob- 
tained. Recently, exponential L 1 convergence to equilibrium was proved in [10] via entropy 
method for global solutions to a simplified convection-diffusion-reaction model with confine- 
ment and Shockley-Read-Hall recombination-generation rate but neglecting the influence of 
the self-consistent potential -0- It would be interesting to study the well-posedness as well as 
long time behavior of the full convection-diffusion-reaction-Poisson system (jl.ip - (jl.2p with 
more general recombination-generation rate R including the higher nonlinear cases (jl.3p and 

(USD. 

For the sake of simplicity we consider the whole-space case posed on M 3 . Similar results 
can be obtained for the two-dimensional whole space case with some minor modifications, 
due to the different properties of the Newtonian potential. Next, we make the following 
assumptions on confining potentials V n ,V p , the recombination- generation rate R and the 
doping profile D: 

(Hla) There exist constants p n ,Pp > such that 

d 2 V d 2 V 

in the sense of positive-defined matrix. Moreover, there exists K > such that 

II a K||l°c(r3) < K, \\AVp\\ L <x>( R 3) < K. (1.6) 

(Hlb) There exists K' > such that 

\\V n (x) - V p (x)\\ L oo m < K' < +oo. (1.7) 
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(H2a) The recombination-generation rate R = R(n,p,x) is of the form 



R(n,p, x) = F(n,p) (np - 5 2 ii n ji p ) 



where ji. 



= e 



V„(x) 



[i p = e 



p( x ) , ^ is a positive constant (the scaled average intrinsic 



it 



carrier density of the semiconductor), which without of loss of generality, can be rescaled 



(H2b) The scalar function F : R 2 — > R is a Lipschitz continuous function with linear growth, 
namely, there exist constants ci,c 2 > independent of n,p such that 

\F(m,Pi) - F(n 2 ,P2)\ < ci(|ni - n 2 \ + \pi — P2I), Vni,pi, n 2 ,p 2 E R, (1.8) 
\F(n,p)\ < c 2 (l + \n\ + \p\) Vti.pel. 

Moreover, F(n,p) > if n,p > 0. 
(H3) D(x) E L^R 3 ) n L°°(R 3 ). 

Remark 1.1. It easily follows from (Hla) that V n (x) and V p (x) are uniformly convex and can 
be bounded from below by a finite number VJ, E M. Thus, ||/x n ||i«= < e Vb and \\HpWl 00 < e Vb ■ 
Without loss of generality, we normalize the confining potentials that 



Remark 1.2. In the following analysis, we set e = 1 without of loss of generality. We note 
that the quasi-neutral limit (namely, zero-Debye-length limit e — > 0) of drift- diffusion- Poisson 
system is a challenging and physically complex modeling problem for bipolar kinetic models of 
semiconductors, which has been analyzed by many authors, see, e.g. [16,23] and the references 
cited therein. 

Now we state the the main results of this paper. 

Theorem 1.1 ( Well-posedness). Suppose that (Hla)-(H3) are satisfied. For any initial data 
m E L 2 (R 3 ,e v '' l ( :r )dx)nL 00 (R 3 ), Pl E L 2 (R 3 , e v " {x) dx) n L°°(R 3 ), m,pi > 0, problem flUJ)- 
(jl.2|) admits a unique global weak solution (n,p,ip) such that for any T > 0, 



n E L oo (0,T;L 2 (R 3 ,e V '" (2:) dx)nL oo (R 3 )), Vn E L 2 (0, T; L 2 (R 3 , e V "^dx)), 
p E L°°(0, T; L 2 (R 3 , e v ^ x) dx) n L°°(M 3 )), Vp E L 2 (0, T; L 2 (R 3 , e v ^ {x) dx)), 
n t ,p t eL 2 (0,T; (H 1 (R 3 ))'), 
n(t) > 0, p(t) > 0, t E [0,T], a.e. a; E R 3 , 
E L oo (0,T;L oo (R 3 )), Ai/> E L°°(0, T; L°°(R 3 )), 



as (5 = 1. 




Besides, we infer from (jl.7p i/tai t/ie norms 
equivalent. 




4 



where ip = tp(t, x) is given by the Newtonian potential 

h( . ,J_f n(t,y)-p(t,y)-D(y) 
ip[l, X) — a I l l y> 

S3 7m3 \x-y\ 

3 

with S3 = ^fjr being the surface area of the 2D unit ball. 
1 (5) 

Theorem 1.2 (Long-time behavior). Under the assumptions in Theorem we have the 
following uniform-in-time estimate for the global solutions: 

IK0llL°°(0,+oo;L>-(R3)) + ||p(t)||l,«(0,+oo;L'-CR 3 )) < °°' Vr E I 1 " +°°]- ( L9 ) 

Moreover, for every fixed t* > 0, the global shifted solution {n{t + s),p(t + s),ip(t + s)) 
(s G (0, i*)) of problem (|l.ip - (|1.2p converges to the unique steady state (uoo , , -i/^oo ) 
satisfies (|4.3p as t ^ +00 in i/ie following sense: 

n(t + -)^^oo, p(t + -)^Poo in L 1 ((0,f) x M 3 ), 
VV(t + -)^V^oo inl^O,**^ 1 ^)), 
^(t + O-'^oo m L 2 (0,t*;L 6 (M 3 )). 

As we have mentioned before, for a class of recombination-generation rate with at most 
linear growth, existence of global weak solutions to problem (|l.ip - (|1.2p has been obtained 
in [23]. However, argument therein fails to apply in our present case due to the higher-order 
nonlinear reaction term R that includes both the band-to-band and the Auger recombination 
(cf. (H2a)-(H2b)). On the other hand, well-posedness results for drift-diffusion-Poisson 
system with higher-order recombination-generation rate have been proved in the bounded 
domain case (see |7-9| for the case with F being bounded, and [26J for the case that F has 
a linear growth). Since now we are considering the Cauchy problem in the whole space that 
is unbounded and the carriers are confined by unbounded external potentials, the methods 
for the initial boundary value problem cannot be used directly. We need to exploit and 
develop several techniques in the literature to prove the global existence and uniqueness of 
solutions to problem (|1.2j) . In order to overcome the difficulties from the higher-order 

reaction term R, we first introduce a L°° cut-off to the unknowns n,p in R and study an 
approximation problem associated with our original system (jl.ip - (jl.2p . To deal with the 
unbounded confining potentials, we then transform the approximate problem into a new 
form by introducing some new variables with proper weight functions. After obtaining the 
well-posedness of the approximate problem, we try to derive proper uniform estimates based 
on a Stampacchia-type L°° estimation technique (cf. [8]) that enable us to pass to limit and 
show the existence of global weak solutions to the original system (|l.ip - (|1.2p . Finally, we 
get uniform-in-time L r (r £ [l,+oo]) estimates for the global solutions under more general 
assumptions by extending the methods in [10, 24J and show the long-time behavior. 

The remaining part of the paper is organized as follows. In Section 2, we prove the 
well-posedness of an approximate problem and obtain some uniform estimates that are in- 
dependent of the approximate parameter. In Section 3, we prove the existence of global 
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solutions to the original problem (ll.ip -( fl~2j) by passing to the limit and show the uniqueness 
of the solution. In Section 4, we obtain some uniform-in-time estimates of the solutions and 
show that as time tends to infinity the global solutions will converge to a unique steady state. 



2 Well-posedness of the Approximate System 

In this paper, we use C, Ci(i G N) to denote genetic constants that may vary in different 
places (even in the same estimate). Particular dependence of those constants will be explained 
in the text if necessary. i2" m (R 3 ) (m £ N) is used to denote the Sobolev space W m,2 ( 
and || • ||#m(R3) is the corresponding norm. We denote L r (M 3 ) (r > 1) with norm || • \\ L r 
Moreover, we denote the vector space L r (IR 3 ) = (L r (R 3 )) 3 (r > 1). The norm on L r (IR 3 ) is 
denoted by || • Hl^r 3 )- For any Hilbert space H, we denote its subspace 

H + = {/(x) G H | f(x) > 0, a.e. x E M 3 }. 

To overcome the difficulty brought by the higher-order nonlinearity R, we introduce the 
following approximate problem. For any a > 0, consider 



(AP)< 

subject to the initial data 



d t n a = div(Vn CT + n CT V(Vv + V n )) - R(n a ,p a ,x), 
d t pa = div(Vp CT +fvV(-Vv + Vp)) - R(n a ,p a ,x), (2.1) 
-A-i/v = n a -p a - D(x), 



n a (x,t)\ t=Q = m, P<j(x,t)\t=o =Pi- (2.2) 
The approximated recombination-generation rate R in (|2.ip is given by 

n a Pa 



R(n a ,p a ,x) = R 



1 + aria- ' 1 + ap a 



v n <r Pa \ n a P<7 \ , OQ x 

= F \T~ 1 >7~; I T~i 7~i AWp J ■ 2 - 3 ) 

V 1 + vn a 1 + ap a J \ 1 + an a 1 + ap a J 

Lemma 2.1. Under assumptions (H2a)~(H2b), the function R = R(n a ,p a ,x) satisfies the 
following properties 

(i) R has at most a linear growth for any n a ,p a > 0, i.e., 

\R(n a ,p a ,x)\ < C a (a(x) + n a + p a ), (2.4) 

where < a(x) £ L^IR 3 ) n L°°(M 3 ) and a(x) £ L 2 (M 3 , e Vn ^dx) n L\ (M 3 , e v ^ x Ux) . 

(ii) R is locally Lip- continuous in L\ (M 3 , e Vi ^dx] , i = {n,p}, such that for any n£\ vffi 
€ L\ (M 3 , e Vn ( x ^ dx) , p { a \ p^ G L\ (M 3 , e v ^dx) , there holds 

\\R{n£\$\ ■) - R{n^J 2 \ OH^^a^w^) 
< K (\\n { ^ - n {2) \\ L2{R3jeVl ( x ) dx) + \\pW -p^ ] || L 2 (R 3 ie vi ( X)dx) ) , (2.5) 

where the constant K may depend on c\,C2,cr and Vj,. 
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Proof. We observe that for any ip > 0, it holds 



<P 1 
< — - — < - and 



<<p, V(T>0. 



1 + acp a 1 + crip 

Due to this simple fact and assumptions (H2a)-(H2b), we can verify that 



(2.6) 



\R(n a ,p a ,x)\ < c 2 1 + 



+ 



pa 



1 + an a 1 + ap a 



Pa 



1 + crn a 1 + ap a 



< c 2 ( 1 + - J finUp + c 2 [ - + \ ) (n CT +Pa)- 



(7 



(2.7) 



Then we can simply set C a = c% (l + ~ -\- J^j and a(x) = fin^p, which obviously satisfies the 
required conditions by assumption (Hla). 

For any i$\ n? G L\ (M. 3 , e v "^dx) , p { }\ p? e L% (R 3 ,e Vp ^dx), we infer from ([231) 
that 



n 



(i) 



(2) 



1 + an? 1 + crni 2) 



<\n^-nf\ 



Pa 



„(2) 
Pa 



1 + ap^ 1 + crpi 2) 



^Ip^-P^l, (2-8) 



and result, 



< 



Pa 



n 



(2) 



„(2) 
Pa 



1 + mi 1 ' 1 + a P a 1] 1 + ani 2) 1 + ap? 



Pa 



1 + ap? 



n 



(i) 



(2) 



1 + otz^ 1 + an? 



+ 



ri 



(2) 



1 + an? 



Pa 



„(2) 
Pa 



1 + crp^ 1 + trpi 2 } 



< -(|n«-n( 2 )| + |pW-p( 2 )|). 

(T 



(2.9) 



Denote 



Then we get 



(?) 



(?) 

Pa 



1 + an?' 1 + ap? , 



, J = 1,2. 



pW ■) - ii(4 2 ),p( 2 ),-)ll L2(R 3^ Wda;) 

< || (i 7 ! - i ? 2)^n^p|| L 2 (R 3 ie ^Wdx) 



+ 



+ 



n 



(i) 



Pa 



n, 



(2) 



„(2) 
Pa 



1 + crn^ 1 + 1 + an^ 1 + up? 



(F x - F 2 ) 



n 



(2) 



„(2) 
Pa 



1 + (rni 2) 1 + a^ 2) 



h + h + h. i = {n,p} 



Mi 



Idx) 



(2.10) 



For the case i = n, it follows from ([22]), ([23]), ([23]), (H2a)-(H2b) and Remark O that 



h < 



Cl 



n, 



(l) 



n. 



(2) 



1 + crn^ 1 + ani 2) 



+ 



Pa 



(2) 
Pa 



1 + api- 1} 1 + apS 2) 
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< 2cie 2 ^ (||n« - || L 2 (K 3 je v nW(ic) + Ib^ - p£ 2) II^s^m^) 



/ 2 < 



+ 



P. 



(i) 



1 + on),' 1 + apX 



(i) 



(l^-^l + l^-pfl) 



Z- 2 (R 3 ,e v "C*)<fa) 



^ ^ + ^ (ll^ - 4 2) |lL2(R3 ie V„ Wd:c) + -^ 2) |lL2 (R 3 ie V nWd2;) ) , 



h < 



< 



(7 

2ci 



Collecting the above estimates together, we see that ()2.5[) (i = n) holds with 



The case i = p can be treated in the same way. The proof is complete. 
Now we state the main result of this section. 



□ 



Theorem 2.1. Suppose that assumptions (Hla)-(H3) are satisfied. For any a > 0, ni G 
L 2 (R 3 ,e v ^ x Ux) n L 4 (M 3 ), pj G L 2 (M 3 , e^dz) n L 4 (R 3 ), n/,p/ > 0, problem ([2TTD - (|2T2l) 
admits a unique global weak solution (n a ,p a ,ij) a ) such that for any T > 0, 



n 



a G C([0,T];L 2 (R 3 ,e v ^dx)), Vn CT G L 2 ((0,T);L 



2/n3 „V, 



Per G C([0,T];L 2 (E 3 ,e^W(i2 ; )), V P(T G L 2 ((0, T); L 2 (M 3 , e^ds)), 

rv(i)>0, p CT (t)>0, tG[0,T], a.e. x G M 3 . 
= ip a {t,x) is the Newtonian potential with respect to n a (t,x),p a (t,x) given by 

1 f n a (t,y)-p a (t,y)-D(y) 



1pa(t,x) 



S3 



\x - y\ 



-dy, with S3 



3 

27T2 



3^ ' 



r(|) 



The proof of Theorem 12.11 consists of several steps. We introduce the following transfor- 
mation of unknown variables (cf. 



u := n a e 2 



v := p a e 2 . 



(2.11) 



A direct computation shows that u and v satisfy the following transformed approximate 
system 

u t - An + A n (x)u = fi(u,v,ip a ), 



(TAP) < 



v t - Av + = f 2 (u,v,ip a ), 

V n (x) Vpjx) 

—I\ip a = ue 2 — ve 2 — JJ(x), 



(2.12) 



S 



where 

A n (x) = \\W n {x)\ 2 - l -^V n {x) + K, A p (x) = \\VV p {x)\ 2 - ~AV p (x) + K. 
System (|2.12|) is subject to the initial data 

Vn(x) v p( x ) 

u(x,t)\ t= o = nj(x)e 2 :=uj, v(x, t)\ t=0 = pi[x)e 2 -= Vl . (2.13) 

It follows from (jl.6p that A n (x) and A p (x) are bounded from below, i.e. 

K , , , K 



A n (x) > A p {x) > a.e. x G 



p3 



Under the transformation (|2.1ip . the right-hand side of (|2.12p are given by: 

Vn(x) ~ 

fl(u,v,ip a ) = Ku + e 2 (div(n CT VVv) - R{n a ,Pa)) 

1 2 VnM VpW 

= Ku + Vu ■ Viper — -uVtpa ■ W n — u e 2 + uve 2 + D[x)u 
-e 2 R(n a ,p a ,x), 

V p (x) 

f2(u,v,ip a ) = Kv + e 2 (div(-p CT VVv) - R(n a ,p a )) 

1 r, _ W 

= Kv — Vv ■ Vip a + -vVip a ■ Wp — ire 2 + uve 2 — D(x)v 
-e 2 R(n a ,pa,x). 

As a first step, we prove the local well-posedness of the transformed approximate problem 
(l2J2l - (f2TT3]l . 



Proposition 2.1. Suppose that (Hla)-(H3) are satisfied anduj, vj G L+(M 3 ). Then for any 
a > 0, i/zere exzsis > suc/i £/iai problem (|2.12p - (|2.13p admits a unique solution {u,v,il) a ) 
on [0, T a \, which satisfies 

u,v G C([0,T a ];L 2 + (R 3 )), Vu,Vv G L 2 (0, T CT ; L 2 (M 3 )), 

uX7V n ,vX7V p G L 2 (0,T CT ;L 2 (IR 3 )). 

T/ie potential ip a is given by 

Va{t,x) = — 1 1 dy. 

J]R3 \x - y\ 

Proof. We consider the following auxiliary linear problem of the transformed approximated 
problem (|2A2]) - p33]l . such that for any u,v G C([0, T]; L 2 (M 3 )), Vu,Vv G L 2 (0, T; L 2 (R 3 )), 
uVV n ,vVV p G L 2 (0,r ; L 2 



(ATAP) < 



u t - An + Ai(a;)ii = (u,v,ijj a ), 

v t -Av + A p (x)v = f+(u,v:,4> r7 ), (2.14) 
u(ar, t) |t =0 = uj, v(x,t)\t= = vi, 
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where ip a satisfies 

-Ai/j a = ue ~-ve ~ - D(x), (2.15) 
and the nonlinearities /{*" , f£ are given by 

f^(u,v,tp a ) = Ku + + e^ 

with u + := max{0,-u}, v + := max{0,v}. 

Since now the nonlinearity R in the approximate problem (|2.ip satisfies the properties in 
Lemma 12.11 using assumptions (Hla)-(H3) we are able to prove the local well-posedness of 
problem f|2.12p — (|2.13|) by the contraction mapping principle as in [241 Theorem 2.2]. Since 
the proof is the same, we only sketch it here. Denote 

S T = {(u,w)GC([0,T];L 2 (R 3 )xL 2 (l 3 )), 

(Vn,V?))Gl 2 (0,T;L 2 (R 3 )xL 2 (R 3 )) : 

n (0) = m > 0, v(0) = vi > 0, max (\\u\\ 2 L2 {R3) + |M|| 2{R3) J < 2M, 
J (l|V^(t)||2 2(R 3 ) + ||Vv(t)||2 2(R3) ) dt < 2M, 

/ / (\u(t)VV n \ 2 + Kt)VK| 2 ) eft < 2M.) 
Jo Jr 3 ' ' J 

where 

M ■= ll«l|||2(R3) + ||U/|| 2 2 (R 3). 

Then we can prove that there exists a sufficiently small T a > such that the mapping 
Q : (u, v) i y (u, v) defined by (|2.14p maps to itself and is a strict contraction. Hence, the 
contraction principle entails that Q has a unique fixed point in St ct such that Q(u, v) = (u, v). 
Next, due to the special structure of the approximated recombination-generation rate R, 
using the idea in [T3], one can show the nonnegativity of the fixed point (u, v) of Q, provided 
that the initial data ui,vj are nonnegative (cf. [241 Theorem 2.2]). Since u + = u, v + = v 
if u, v > 0, we see that f^(u,v,ip a ) = fi(u,v,ip a ) (i = 1,2) for u, v > 0. Thus, (u,v) is the 
local solution of problem (|2.12p - (|2.13p . The details are omitted here. □ 

Lemma 2.2. Assume that (Hla)-(H3) are satisfied. For any T > 0, if nj,pj E L r (M 3 ) n 
L 1 (M 3 ) ; re N, we have 

IMt)|| L .(tt») + ||p<r(t)||L.(R3)<C , r , 0<t<T, s = l,...,r. (2.16) 

Moreover, 

IIV^Hloo^) <C t , 0<t<T, (2.17) 

provided that nj,pi £ L 4 (R 3 ) n L X (]R 3 ). In particular, the constant Ct is independent of 
a>0. 



div ( u e 



R u + e~ 



Vn(x) 



, v 1 e 



2 , X 



-div 



v e 



R 



u e 



2 



v e 



Vp(g) 
2 , X 



10 



Proof. Integrating the equations for n a and p a in (12.1H on M 3 , we infer from (|2.6|) . (H2a) 
(H2b) that 



(it 



n a Pa 



< 2 / /i^/Xp-F 



(1 + an a )(l + crp CT ) 

"(j Per 



1 + an<j ' 1 + ap a 

< 2c 2 / /U n /i p (l + n CT +p fT )dx 

< C(l + ||n CT || L i (M 3) + Hpo-Hh 



dx 



Pa 



1 + (771(7 1 + ap a 
n a Pa 



1 + an a 1 + ap 



, x ) dx 
dx 



(2.18) 



which yields 

IK(*)IIl1(R3) + l|2V(£)|Ul(R3) < (||^/||l1(K 3 ) + IIP/Ill, 1 ^ 3 ) + 1) 



ViG[0,T]. (2.19) 



Next, multiplying the equations for n a and p a by n r a ,p r a (r £ N) respectively, integrating on 
M 3 and adding the resultants, we obtain that 



1 d 
r + ldt 



n r + 1 +p r + 1 )dx + 



4r 



\r + lf 



r+l 



r + l 
2 



da; 



r + 



(2.20) 



i? 



P<7 



1 + (TTLq- 1 + crpo- 
Since n a ,p a > 0, we observe that 

™<r Pa 



,x) (n r a +pl)dx. 



R 



1 + an a ' 1 + <?Pa 
n a Pc r- 



X ) ( n a+Pa)dx 
Pa 



(1 + ana)(l + apa) \l + ana 1 + apa 

+ / p n p p F ( — — _EZ ) („r +^)rfx 

7 R 3 ^ V 1 + an a 1 + VPaJ 



« +Pa)dx 



(2.21) 



<C f {nl+pl + rf+ l +p r + l )dx. 
On the other hand, from the Poisson equation and the elementary calculation 

r 

( Q r+i _ 6 r+i^ 6 _ a) = _ ^ a r ~ l b k (a - bf < 0, Va, b > 0, 



(2.22) 



fc=0 



we infer from (H3) that 

-J- / A^K +1 -^ +1 )dx 

^ + 1 JR3 



r + l 



(p (T -n <T + Z?(x))« +1 -^ +1 )^ 



< p(x)|| ioo(K3) / (n^+p^dx. 



(2.23) 
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Besides, it follows from (Hla) that 

r - [ (AV n n r + 1 + AV pP r a +1 )dx < K I « +1 + p r + 1 )dx. (2.24) 



r + 

Summing up, we have 
1 d 



r + ldt 
C 



n r a +1 +p r + 1 )dx + 



4r 



{r + lf 



r + l 

V(n CT 2 



+ 



r+l 



dx 



n 



r + pl + nl +1 +p r + 1 )dx 



< C I {l + n r a +1 +pl +l )dx. 

Then it follows from the Gronwall inequality that 

cT 



(2.25) 



n 



I r+l 
a\\ jjr+i 



+ l|Po-|l2t+i(K3) + 



JUL 3 



r+l 

V(n 2 



+ 



V [p 



r+l 
2 



dxdt <C T , r G N, 
(2.26) 



provided that J R3 (n V j +1 + p T j +1 ) dx < oo. 
Since 

1 f n a (t,y) -p a (t,y) - D(y) 



<->3 



dy, 



we have 



For any x G 



Vijj a (t,x) 



S 3 



|V^(t,x)| <C 



F - 2/| 

(ry(t, y) - y) - D{y))(x - y) 
\x — y\ 3 

\n<r(t,y) -Pa(t,y) - D(y)\ 



dy. 



\x - y\ 



dy. 



(2.27) 
(2.28) 

(2.29) 



It follows from [191 Corollary 2.2] (a direct consequence of the Hardy-Littlewood-Sobolev 
inequality [27j ) that for any q, q' G (l,+oo) with - = — ^, 

||VV>er||L9(K3) < C||n CT - Z)(x)|| Xa / (R3 )- (2.30) 

Besides, if n^p^ G L°°(0, T; L 9 '(M 3 ) n L X (]R 3 )) with > 3, then we have 

Vie^^TiL^tR 3 )). (2.31) 

The proof is complete. □ 

Based on the a priori estimates obtained in Lemma 12.21 we can prove existence of global 
solutions to problem (prL|) - ([2T2"|) . 



Proposition 2.2. Suppose that all assumptions in Proposition \2.1\ are satisfied. Assume in 
addition that nj,pi G L 4 (1R 3 ). The local solution (n a ,p a ,ip a ) obtained in Proposition \2.1\ is 
global. 
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Proof. Multiplying the first two equations in (|2.12)) by u and v respectively, integrating on 



we get 



1 d 112 



+ ||V«| 



L 2 (R3) 



+ 



I M 



x)u 2 dx 



fi(u,v,i() a )udx 



K I u 2 dx + - I u 2 Aip a dx - i / u 2 Vif> a ■ VV n dx 

2 ./ra>3 2 



- / ue 2 R 



ue 2 



ve 2 



1 + erne 2 i _|- ave 2 



, x dx, 



(2.32) 



li 2(R 3) + || Vw||l2 {R 3) + J A p (x)v 2 dx = J f 2 (u, v, ip a )udx 



Id. 
2dt' 

/v / v 2 dx - - / v 2 A^ a dx + - / v 2 V^ CT • W p dx 
2 7to3 2 



ve 2 R 



ue 2 



Vp(s) 

2 



, 1 + aue 2 i _|_ ave 2 



(2.33) 



From the Poisson equation for ip a , Lemma[2]2]and the Gagliado-Nirenberg inequality ||n|| s 



< 



C||V«| 



L 2 (R3)II U Hl 2 (R3)' We haVe 



< 



< 



/ (u 2 - v 2 )Aip a dx = / (u 2 - v 2 )(-n a + p a )dx + / D(x)(u 2 - v 2 )dx 

P(^)lk~(R3)(lMl!,2( R 3) + IhH^^S)) 



1 



v „ ■ , , + l|Vv||^ (R3) J +C[\\u\\{ 2m + |b||i 2(K3) J. (2.34) 
We infer from the assumption uj,vj E L^R 3 ) and the Cauchy-Schwarz inequality that 



1 



1 



1 



1 



njdx < — / /i n dx + — / u^dx < +oo, / p/dx < — / « p dx + - / f/dx < +oo, 

2 JR3 2 J R 3 J R 3 2 J R 3 2 J R 3 

(2.35) 

namely, nj,pi E L 1 (M 3 ). Lemma 12.21 implies that if the initial data satisfy ni,pi E _L 4 (M 3 ) Pi 
L X (IR 3 ), then the estimate (|2. 17|) holds. As a consequence, we get 



u 2 VVv • W n dx 



< l|V'0 o -|| LoO ( M 3)||'u|| I/ 2 (R 3)( / u 2 |W n | 2 dz 

<\ I « 2 |VK| 2 dx + C T ||n 
2 Jr3 



(2.36) 



and 



v 2 Vip a ■ Wpdx 



< IIV^Hl^Ors)!^ IU 2 (r 3 )( / v 2 \W p \ 2 dx^ ' 

1 [ W 2 |W p | 2 dx + C T |b 

2 JR3 



(2.37) 



< 
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It follows from the nonnegativity of (u,v) and assumptions (Hlb), (H2b) that 

x / YaM _VpW \ 

ue 2 +ve 2 Ji? ZW ZW X ) 

' \1 + aue 2 1 + ave 2 / 

/ / j^W VpM\ / 2 2 \ 

< / ue 2 + ve a I F I ™, ^yinnfxpdx 

V / yi + o-ue 2— l + ave ~ / 

_ f ( YbM V p(x) \ ( VnW _W\ 

< C J I ue 2 -\- ve 2 \ I u n u p + ue 2 + we 2 \ dx 

< C^ + Hl^aj + IHli^a)). (2.38) 
As a result, we have 



+ / u 2 \VV n \ 2 dx + [ v 2 \VV p \ 2 dx 
< C T (l + \\u\\ 2 L2m + \\v\\ 2 L2m ). (2.39) 
It follows from the Gronwall inequality that for any T > and t G [0, T], 

II«(*)IU»(rs) + IK*)lb( R 3) < Or, (2.40) 
/ (UVull^^s) + ||V^||^ (R3) )dt < C T , (2.41) 

/ / (u 2 \VV n \ 2 dx + v 2 \VV p \ 2 )dxdt < C T , (2.42) 
Jo Jr 3 

where Op is a constant depending on T,c 2 , V b , ||u/|| L 2 {R 3), ||v/||x,2(M3), ||"/||l 4 (r3), IIp/||l 4 (r3) 
but it is independent of a. 

Then we are able to extend the local solution (n a ,p a ,ip CT ) obtained in Proposition 12.11 to 
the interval [0, T] for arbitrary T > 0. The proof is complete. □ 

Proof Theorem 12. 1L Recalling the transformation (12. lip , we obtain the conclusion in 
Theorem 12.11 from Propositions 12. 1| 12.21 

3 Well-posedness of the Original Problem 

In this section, we prove the existence and uniqueness of global solutions to the Cauchy 
problem of original system (jl.ip . For this purpose, we shall derive some a priori estimates 
on the solutions (n a ,p a ,i/j a ) of the approximate problem that are uniform in the parameter 
a > 0. Then, we pass to the limit as a — > to achieve our goal. In Lemma [2.2l we have already 
shown the uniform estimates on ||w<r(i)IUf(K3)) l|Po-(0llL r (R 3 ) 011 arbitrary interval [0, T]. Next, 
we proceed to obtain uniform estimates on the L°° norms of n a and p a via a Stampacchia- 
type L°° estimation technique (cf. [8]). The following technical lemma (see [H Lemma 2.3]) 
plays an important role in the proof. 
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Lemma 3.1. Suppose u(k) is a nonnegative non-increasing function on [Aft, +00), and there 
are constants 7, /3 such that 

w(jfe) < M(k)(k- fc)- 7 u;(fc) 1+/3 , \/k>k>k , (3.1) 

where the function M{k) is non- decreasing and satisfies 

< k" 1 M{k) < M , V e [fco, +00). (3.2) 

T/ien 

w(fc*) = 0, u*ifc fc* = 2A ll + 2^M ' 37 w(fc )~ ) • (3.3) 

Lemma 3.2. Suppose that all assumptions in Theorem \2.1\ are satisfied. Assume in addition 
that ni,pi £ L°°(M 3 ). Then for any T > 0, we have 



\\n a (t)\\ Laom <C T , ||Pa(t)IU»CR3)<Qr, 0<t<T, (3.4) 
where the constant Ct is independent of a > 0. 
Proof. Denote 

k : = max{||n/|| LO o (R 3), ||p/||l-(r3)} > 0. (3.5) 

For any k > ko, we set 

B nk (t) = {x£R 3 ,n a (t,x) >k}, B pk (t) = {x£R 3 , Pa (t,x)>k}. (3.6) 
and for arbitrary T > 0, 

ur(k) = sup (\B nk {t)\ + \B pk {t)\). (3.7) 

0<t<T 

It is obvious that wr(fc) is a nonnegative, non-increasing function on [feo,+oo). Moreover, 
the uniform bound (|2.19p implies that u>t(A) is bounded for any T > 0. 

For any k > ko, multiplying the first and second equation in ()2. 1|) by (n a — k) + and 

(p a — k) + , respectively, integrating on M 3 and adding the resultants together, we obtain 



1 d 



(lIK - k) + \\ 2 L 2 m + \\( P a - /c) + ||| 2(R 3)) + ||V(n CT - k) + \\l 2 



2dt 

II 2 



+\\v(p a -k) + " 2 



= \ [dw(n a V{ip a + V n ))(n a -k) + + div(p a V(-ip a + V p ))(p a -k) + ]dx 
Jr 3 

~[ R( T ^^, T ^^^)[(n a -k) + + ( Pa -k) + ]dx 

J R 3 V 1 + an v 1 + G Vcj ) 

=: h+h- (3.8) 

Integrating by parts and using the Poisson equation for ip a , we expand the first term I\ as 
follows 

h = k f [A(^ + V n ))(n a -k) + + A(-^ + V p ))( Pa -k) + }dx 
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+ / (div[(n ff - fc)+V(Va + V n )](n a - k) + + div[( P(J - k) + V(-^ a + y p )](p ff - k)+) dx 

fc / + ^ + [(n ff - k) + - (p a - k) + ] dx 

it 3 

+A; / (AV n (n a - k) + + A^( P(7 - k)+) dx 

JR 3 



v(w + v n ) ■ v[K - k) + ydx - - / v(-v CT + v p ) ■ v[(p ff - AO+rdx 

2 JR3 

/ (— n ff + [(?V - k) + - (p a - k) + ] dx 
Jr 3 

+k / D(x) [(n a - k) + - {p a - k) + ) dx 
Jr 3 

+k [ (AV n (n a - k) + + AV p (p a - k)+) dx 
Jr 3 

+ 1 I {-n a +p a )([(n a -k) + ] 2 -[{p a -k)+] 2 )dx 
1 Jr 3 

+~ I D(x)([(n a - k) + ] 2 - [(p a - k) + ] 2 )dx 

* Jr 3 

+1 [ (AK[(n CT - k) + ] 2 + AV p [( Pa - k) + ] 2 )dx 

* Jr. 3 



=: Ji + ... + J fl . 
It is easy to verify that 



(3.9) 
(3.10) 



Ji < 0, J 4 < 0, Mk > k . 
Besides, it follows from (H2b) that 

h < C2 / MnMp(l + n a +Pa)[(n a - k) + + (p a - k) + ]dx := J 7 . (3.11) 

Jr 3 

Integrating (|3.8p with respect to time, we infer from (|3.9p . (|3.10p and (|3.1ip that 

\ (||K(t) - k) + \\ 2 L2m + \\{p a (t) - k) + \\ 2 L2m ) 



+ / (HVK-A;)- 



+ \\V(p a -k) + \\l 2(m )dT 



< [ (J2 + J 3 + Js + Je + Ji)dr. 
Jo 



(3.12) 



By assumptions (Hlb) and (H3), we have 



J^dr 



+ 



J§dr 
t 



< -\\D(x)\\ Lxm J \\(n a -k) + \\ z L2im + \\(p a -k) + \\ z L2{R3) dT 
+ \\AV n \\ Loo{R3) \\{n a - k) + \\ 2 L2{R3) dT + \\AV p \\ L o 0(R 3 ) ||( 
\(n a - k) + \\ 2 L2(R3) + \\(p a - k) + \\ 2 L2(R3) dT. 



< C 



k ) + \\h(R 3 ) dT 

(3.13) 
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Furthermore, using the Holder inequality and Gagliardo-Nirenberg inequality, we get 



J 2 dr 



+ 



J 3 dr 



< k f I (\D(x)\ + \AV n \)(n a -k) + dxdT + k f f (\D{x)\ + \AV p \)( Pa -k) + dxdT 
Jo Jr: 3 Jo Jr 3 



< k(\\D{x)\\ Loo0IL3) + \\AV n \\ Lo 



+k(\\D(x)\\ L ^ m + HA^IIioo^s)) 



<x - k) + \\ L 3( R 3)\B nk \*dT 

, 2 

Pa ~ k) + ||L3(]R3)|Bp fe |3dr 



dr 



2 r i i 

< Cku T (k)3 \\V(n a - k) + \\l 2m \\(n a - k)+\\ 2 L2m 

+Cku T {k)l j* ||V(p ff - k) + \\l 2m \\{ Pa - k) + \\l m dr 

< V T j* ||VK - k) + \\ 2 LHm + ||V(p ff - A;) + ||2 2(K 3 ) dT 
+r/T^ UK - k) + f L2m + UK - k) + \\ 2 L2m dr 



Next, 



C 4 

H — k 2 oj T (k)3. 

7] 



J-jdr 

o 
t 



(3.14) 



< [ [ ^ P [{n a -k) + + {p a -k) + +2k + l}[(n a -k) + + { Pa -k) + ]dxdT 
Jo Jr 3 

< Cj* ||K " ^) + lli 2 ( R 3 ) + UK " k^Wl^dr 

+C{2k + l) / / K - k) + dxdr + C(2k + 1) f [ { Pa -k) + dxdT 
Jo Jr 3 Jo Jr 3 

< cj o (||K - k) + \\h m + ||K - k) + \\ 2 L2m )dT 

+ C(2k + l)u T {k) 2 S [ (||K " fe) + IU3( K 3) + ||K - k) + IU3( R 3)) dT 

Jo 

< rfrj \\V(n a -k) + f L2m + \\V(p a -k) + \\l 2m dr 
+(C + tfH J* ||K - fc) + ll! 2 (R3) + ||K - k) + \\ 2 LHR3) dr 



c 

+-(2k + l) 2 uj T {ky 



(3.15) 



17 



Taking r\ = in the above estimates, we obtain that 

\\(n a (t)-k)^l 2m + \\(p a (t)-kni 2m 
+ [ (\\V(n a -k) + \\ 2 L2m + \\V( J)a -k) + \\l 2m )d7 



o 



(3.16) 



+ l|2 



dr 



<Ci J (\\Mr)-k)+\\l 2iR3) + \\{p a (r)-k) 

+ TC 2 (k 2 + l)ur( k) t. 

It follows from the Gronwall inequality that 

- fc)+||2 2(R3) + \\( Pa (t) - k) + \\ 2 L2m < e c ^ T TC 2 {k 2 + l)co T (k)l, V * G [0,T]. 

(3.17) 

On the other hand, for any i £ [0, T] and k > k > ko, 

\\(n a (t)-k) + \\ 2 L2(R3) = [ \(n a (t,-)-k) + \ 2 dx> [ (n a (t,-)-k) 2 dx 

JB nk (t) J B nk (t) 

> (k - k) 2 \B nk (t)\, (3.18) 

\\(pa(t)-k)+\\ 2 L2m = f \( P a(t,-)-k) + \ 2 dx> f ( Pa {t,-)-kfdx 

JB pk (t) JB pk {t) 

> (k-k) 2 \B pk (t)\. (3.19) 
We deduce from (pTT7]) - (prT9"]) that 

uj T {k) <e ClT TC 2 {k 2 + \){k-k)- 2 uj T {k)^, Vk > k > ko- (3.20) 
Now in Lemma |3. 11 we set 

M(k) = e ClT TC 2 (k 2 + 1) > 0, 7 = 2, (3 = \. (3.21) 
The function M(k) has the following property 

^ = ^e c ^TC 2 < (l + ^) e ClT TC 2 := M , V fce ^,+oc). (3.22) 
Therefore, there exists a constant 

>15 



= 2k (1 + 2 iO M ^ T (/c )3j > ko, (3.23) 
which is independent of a such that 

w r (£;*) = 0. (3.24) 

Namely, 

n a (t,x)<k* and p a (t, x) < k* , a.e. in M 3 x [0, T]. (3.25) 

The proof is complete. □ 
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Proof of Theorem ll.lL Lemma f2.2l and Lemma f3.2l yield the following uniform estimates 
that are independent of the parameter a > 0: 



\ n a\\L°°(0,T;Li(R 3 )) + \\Pcr || L°° (0,T;L«(R 3 )) < Or, 
II ? V||l 2 (0,T;/P(IR 3 )) + 11^1^2(0^.^1(^3)) < Ct, 



(3.26) 
(3.27) 



l|V^<7|| L oo( 0iT . L ,'( R 3)) + || A^ CT IIloo^l^irS)) < C T , <?e(0,+oo), (3.28) 
where q £ [l,+oo], g' G (l,+oo]. Besides, we infer from (12.42p . (12. lip and (HI) that 

I [ (nl\VV n \ 2 dx+pl\VV p \ 2 )dxdt <e~ Vb [ [ (u 2 \VV n \ 2 dx + v 2 \ VV p \ 2 )dxdt < C T . 
Jo Jm. 3 ' Jo Jm 3 

(3.29) 

Then by the equations for n a and p a in (|2.ip and (|3,26p - (l3.29p . we obtain that 



lldiTVllL^O.T;^ 1 ^))/) + ||%J CT ||^2 (0,T;(.ffl (R 3 ))') < ^Tj 

From the uniform estimates f|3.26p — (|3.28|) . we deduce that there exist 

n,p£ L 00 ^,^ 00 ^ 3 )), 

with 



(3.30) 



V M ,Vpei 2 (0,r ; L 2 (I 3 )), ^^GL^OJ^ff^R 3 ))'), 

and V with A^ G L°°(0, T; L°°(IR 3 )), G L 2 (0, T; L 2 (M 3 )) such that for a sequence {cr,-} 
as j — > +oo (not relabeled when taking a subsequence), 



(3.31) 
(3.32) 



n a . ->> n, p aj 



p, weakly-star in L°°(0, T; L°°(M 3 )), 
Vn CTj Vn, Vp CTj . Vp, weakly in L 2 (0,T; L 2 (M 3 )), 
d t n a . d t n, d tPaj -» %>, weakly in L 2 (0, T; (^(M 3 ))'), 

A^, weakly-star in L°°(0, T; L°°(IR 3 )), 
V^ CT . ->> V^, weakly in L 2 (0, T; L 2 (M 3 )). 

Moreover, on account of the compact embedding theorem we have 

p, strongly in L 2 (0,T; L 2 oc (R 3 )), thus also a.e. in (0,T) x 
00^3^ we have 



(3.33) 
(3.34) 
(3.35) 
(3.36) 
(3.37) 

(3.38) 



Then, for any (p <G L 2 (0, T; C, 



< 



Jm 3 
T 



V^ ffj — nVip)Vipdxdt 



/ / (n CT - n)Vi) a N^pdxdt 
Jo Jm 3 



+ 



(Vt/v, - V*p)nV<pdxdt 



< || VV'ctj Hl^CTjL^R 3 ))!!™^ - HIl 2 (0,T;L 2 (sup PV ))II V ^IIl 2 (0,T;L°° 

+ [ [ (Vipa- -Vip)(riVip)dxdt 
Jo Jm. 3 

— > 0, as Oj — >• 0, 



(3.39) 
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— > 0, as Oj — > 0. 



(3.40) 



similarly, 

Jo Jrz 

Next, we study the convergence of the recombination-generation rate. The uniform bound 
([336]) and (H2b) yield that 



< C T . 



L 2 (0,T;L 2 ( 

Thus there exists G £ L 2 (0, T; L 2 (R 3 )) such that (up to a subsequence) 

R(n a , PcT ,x) -> G, weakly in L 2 (0, T; L 2 (R 3 )) as Oj -> 0. 



(3.41) 



(3.42) 



For any bounded domain f2 C M , there holds 
' ■ 1 ' n n 



o Jn 



1 + Ojn (Tj 
which implies that 



\ 2 




J dxdt < 


fj 




Jo Jtt 



(n a — n) 2 dxdt + 



o Jn 



(ojn aj n) 2 dxdt 



0, as Oj — > 0, 



1 + OjU (Ti 



n strongly in L 2 (0, T; L 2 oc (R 3 )) 



In the same manner, we have 

Pa, 



1 + CTjPaj 



p strongly in L 2 (0, T; L 2 oc (R 3 )). 



(3.43) 



(3.44) 



(3.45) 



Since F is Lip-continuous (see (H2b)), we infer from (|3,44l) and (13.450 that on any bounded 
domain Qcl 3 , 



o Jn 

T 



Po 



1 + T;//„. 1 + OjP aj 

2 



< c 



n n 



1 2 



F(n,p) 



o Jn \l + (Tjn aj 



n + 



Pa, 



1 + OjP aj 



dxdt 



p I dxdt 



0, as Oj — > 0, 



(3.46) 



namely, 



F 



Pa, 



F(n,p) strongly in L 2 (0,T;L 2 OC (IR 3 )) and a.e. in (0,T)xR 3 . 

(3.47) 



(3.48) 



l + CTjre CTj l + Ojp a 

As a result, we have the point-wise convergence of R 

R(n a ,Pa,x) — > R{n,p,x), a.e. in (0,T) x M 3 , 
which together with (|3.42p implies that G = R(n,p,x) and 

R(n a ,p a ,x) ^ R(n,p,x), weakly in L 2 (0, T; L 2 (M 3 )) as Oj 0. (3.49) 
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Based on the above convergent results, now we are able to pass to the limit by letting Oj — > 
in the approximate problem (j2.ip and obtain a global solution (n,p, ip) of problem (jl.ip - (jl.2p . 
The system (jl.ip is satisfied in the following sense that for any <p £ L 2 (0, T; H l ( 



j (n t ,ip) u-i jjidt + / / (Vn + nV{ip + V n ))Vipdxdt + / / R(n,p,x)ipdxdt = 0, 
Jo Jo Jr 3 Jo Jr 3 

(3.50) 

/ (Pt,<p)H- 1 ,mdt+ / (Vp+pV(-ip + V p ))V(pdxdt+ / R(n,p,x)(pdxdt = 0, 
Jo Jo Jr 3 Jo Jr 3 

(3.51) 



VV> • Vipdxdt 



(n — p — D)ipdxdt. 



(3.52) 



Finally, we prove the uniqueness of global solutions to problem (ll,ip - (ll.2p . Let (rn,pi, V>i) 
(i = 1,2) be two solutions to problem (jl.ip - (jl.2p with initial data nn,pn. Set n = ri\ — ra 2 , 
p = p\ — p2, tp = ipi — tp2, tii = nn — ni2 and pj = pn — pi2- Taking the difference of the 
equations for ri\ and 712, testing the resultant by n, we find that 

ft 



1„ 

— n 

2" 

1„ 



+ 



iVn 



-.dr 



L 2 (R3) C 
nVn • WndxdT 



- \R(ni,pi,x) - R(n 2 ,P2,x) )ndxdr 

Jo Jr 3 v ' 

1 2 

Using the uniform estimates for rii,pi,ipi, we have 

|-El+-E2| < ^||AK||io°(K3) + ||A^l||jr,cx.( R 3)^ J 



(riVipi + H2^/ip) Vndxdr 



(3.53) 



dr 



+ ll"2||i°°(0,T;L 3 

Then by (H2a)-(H2b), we deduce that 



e||Vn||£ 2(R3) + C £ ||VV||£6 (M 3) j dr. (3.54) 



Jo 


JR 3 


+ 


Jo Jr 3 


+ 


ff 

Jo Jr 3 


c( 


V b , \\rii 


X J 


An 





(F(m,pi) - F(n2,p 2 ))nfi n fi p dxdT 
(F(ni,pi) - F(n2,P2))nipindxdT 
F(n2,P2){iT'ip + np2)ndxdT 



\L°°(0,T;L°° 



) 



+ 



)dr. 



(3.55) 



In a similar manner, we have the following estimate for p 



+ 



\\Vp\\ 







dr 
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< \\\pi\\l 2 (r3) + - [\\AV p \\ Lx[r3) + \\A^i\\ Loa 



\p\\ 2 L 2 (R3)dr 



(It 



+ INIl°°(0,T;Z,3(R3)) J ( e ll V PllL2(R3) +CyVV>|| L 6 (R 3 ) 
+C (Vb, ||n»||i«>(0,r;L«>(Ill3))j lbi|li°°(0,T;L°°(M3))) y (|M||2( R 3) + ||p|||2( R 3))dr(3.56) 

Since ip satisfies the Poisson equation — = n — p, then it follows from |19t Corollary 2.2] 
that 



Therefore, taking e sufficiently small satisfying 



||W>||l6 (R 3) < C (\\n\\ 2 L 2 (m + ||p||| 2 



.Hp-"-' 



< e < — min •! 1, \\ri2 



1-1 



lL°°(0,T;L 3 



2 IIl°°(0,T;L 3 



}■ 



(3.57) 



(3.58) 



we deduce from (|3.53p - (|3.57p that 



\n(t)\\ 



L 2 



+ 



L 2 



+ / i II V/'/ 1| 7,, i, + \\Vp\\l. 

2 . n„ 1 1 2 , n II 11 11 2 , 11 11 2 



(It 



^ IKIIl 2 (R3) + \\pi\\ L 2 (R3) +C T I [\\n\\ L 2 (R:)) + \\p\\ L 2 



(3.59) 



dr. 



From the Gronwall inequality, we can conclude the continuous dependence on the initial data 
that 
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I-I.v 1 ) T H/' M Vlli2( R 3) < (ll«/|ll2( R 3) + Ib/||12( R 3)J e° Tt , V* G [0,T], 

which yields the uniqueness. The proof is complete. 



(3.60) 



4 Remark on the Long-time Behavior 

In the last section, we have proved the existence and uniqueness of global solutions to 
problem (|1.1|) - (|1.2|) . However, the global-in-time estimates for the solution (n,p,ip) depends 
on T that can be chosen arbitrary. In this section, we extend the results in |1CHI24| to our 
current problem f)l.l|) — ()1.2|) . For this purpose, we first need to obtain some uniform-in-time 
estimates on the global solution. 

Let a = J R 3(nj — pj)dx. We easily see from (jl.ip that the difference of mass is conserved 
for all t > 0: 

(n(t, •) -p(t, -))dx = a. (4.1) 



/ 



The relative entropy associated with (jl.ip is as follows: 



e(t) 
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n In 1 + n 



dx + 



p ( In — - 1 ) + poo 

Poo 



dx 



+ 1 I IW - V^oo\ 2 dx, 

2 JR3 



(4.2) 
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where (noo,£>ooj V'oo) is the steady state of system (jl.lh that satisfies 



n 00 (x) = D n e~^°° u n , A, G R+ 
PoD (x) = Dpe^fip, D n e R+, 

l^ooPoo — l^nl^pi Jjj3 Tloadx Jjg>3 p^dx — (X, 
-Alpoo = Hoc - Poo - D(x). 



(4.3) 



Remark 4.1. Denote I = f R3 e Vn ^dx,J = f R3 e^°° Vp ^dx. Then the coefficients D Ti 
and Dp in fj4.3j) are given by (cf. \2J\ Lemma 3.1]) 



Dr. 



a 



+ Va 2 + 4JJ 



Va 2 + 4JJ 



v ^ -a 
21 ' p 2J v ; 

Following the argument in [241 Theorem 3.1], where the special case V n = V p = V was 
considered, we can still prove the existence and uniqueness of (riooiPoo^oo)- 

Proposition 4.1. Suppose that assumptions (HI a), (Hlb) and (H3) are satisfied. Then the 
stationary problem (|4.3p admits a unique solution (V'oo; ^007 Poo) such that ipoo & D ' 1 
{(j){x) E L 6 (R 3 ) I E L 2 (M 3 )}. Moreover, ^ E L°°(M 3 ) and VV>oo G L°°(K 3 ) n L 2 ( 

The we have 

Lemma 4.1. Suppose the assumptions of Theorem \1.1\ are satisfied. The global solution 
(n,p) of problem (ll.ll) - (ll.2p satisfies 



sup M|n 



+ 



< 00, Vr E [1, +00]. 



(4.5) 



Proof. By a straightforward calculation, we have the dissipation of the relative entropy 



j t e{t) = -V(t) < 0, 



with the entropy dissipation 



D(t) 



/ 



n 



Vlnl- 



n 



dx 



Vln[£ 



R(n,p, x) In 



where N = D n e~ m P = D p e m 



(4.6) 



dx, (4.7) 



(4- 



Based on the entropy dissipation inequality, we can obtain uniform Z7 bounds (r E [1, +00)) 
for n(t) and p(t) exactly as in [24^ Lemma 4.1]. It only remains to show the uniform L°° 
estimate. We note that L°° bounds of solutions to a simplified drift-diffusion system (without 
self-consistent potential ip and with a recombination-generation rate of Shockley-Read-Hall 
type) have been obtained in |10| via a Nash-Moser type iteration method and the results could 
be extended to the case with self-consistent potential [11]. For the convenience of the readers, 
we sketch the proof for our present case with a more general recombination-generation rate. 
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For r > 2, using integration by parts and the nonnegativity of n,p, we get 



d_ 
~dt 



n r+1 + p r+1 )dx + 



4r 
r + 1 



VI n 2 



+ 



V p 



r+l 



dx 



A^{n r+1 - p r+l ) 



+ r [ (AV n n r+1 + AV p p r+1 )dx 
— (r+l) / R(n,p,x)(n r + p r )dx 
= -r [ (n - p){n r+1 - p r+1 )dx + r [ D(x){n r+1 - p r+1 )dx 

JR 3 JR3 

+r f (AV n n r+1 + AV p p r+1 )dx 

-(r + l) / F(n,p)(n r+1 p + np r+1 )dx+ (r + 1) / iJ 2 F(n,p)(n r + p r )dx 

< KII-OIUoofRS) + II a K||l°°(r3)) / nr+1 ^ + r '(ll- D llL-»(R3) + ||AV^|| LO o (]R 3)) / p r+1 cix 

JR 3 JR 3 

+C(r + 1)/ (1 +n + p)(n r +p r ')<ix 
./r 3 

< C(r + 1) / (n r+1 + p r+1 )cfe + C(r + 1) / (n r + p r )dx, 

JR 3 JR3 



where we use the facts that 



Since 



{n-p){n r+l -p r+1 ) > 0, F(n,p){n r+1 p + np r+1 ) > 0. 



2 — r 2— r 

(r — l)r r - 1 > for r > 2, lim (r — l)r'— 1 = 1, 

r— >+oo 



(4.9) 

(4.10) 
(4.11) 



by the Young's inequality and the uniform L 1 estimate of n,p, we deduce that 

/ (n r +p r )dx < 4t / (n + p)dx + (r - l)r^ / (n r+1 + p r+1 )(ix 
Jr 3 ^ Jr 3 Jr 3 

< f (n r+1 +p r+1 )dx, Vr>2, 

^ 7r 3 



(4.12) 



where the constants C, C" are independent of r. Therefore, it follows from (|4.9p and ()4.12j) 
that 



(n r+1 + p r+1 )dx + 



4r 
r + l 



V(n 2 



+ 



V (p 



r + l 



dx 



< C{r + 1) f (n r+1 + p r+l )dx + -, 
7r 3 r 



(4.13) 



where C is independent of r. Based on the differential inequality (I4.13h . we can argue as 
in [10\ Supplement Lemma 5.1] to obtain the uniform L°° bounds for n(t) and pit). The 
proof is complete. □ 



Proof of Theorem [TT2l Lemma [4.11 yields the uniform-in-time L r estimates (|1 .9|) . Then 
the conclusion of Theorem 11.21 follows from the same argument as in [24^ Theorem 4.1]. The 
proof is complete. 
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